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We study the electric-field induced capillary interaction 
of charged particles at a polar interface. The algebraic tails 
of the electrostatic pressure of each charge results in a de- 
formation of the interface u ~ p~ 4 . The resulting capillary 
interaction is repulsive and varies as p~ 6 with the particle 
distance. As a consequence, electric-field induced capillary 
forces cannot be at the origin of the secondary minimum ob- 
served recently for charged PMMA particles at on oil-water 
interface. 

Electrostatic forces operating on charged particles de- 
termine to a large extent the properties of emulsions and 
foams, and stabilize colloidal aggregates [1,2]. Charged 
particles at interfaces or in thin films form 2D hexagonal 
crystals [3-7] , while for bulk colloidal suspensions several 
3D crystal phases have been reported [8]. At a typical 
lattice spacing of a few microns, the interparticle forces 
are of the order of picoNewton. 

A charge at a polar interface and its counterion cloud 
carry a finite dipole moment perpendicular to the inter- 
face (Cf. Fig. 1.) The resulting electrostatic force act- 
ing on two neighboring particles is repulsive and varies 
with the inverse fourth power of their distance [3,9-12]. 
This law has been confirmed experimentally for charged 
polystyrene particles [6] . 

The pair potential of PMMA particles at an oil-water 
interface has been measured by recording the trajectories 
and evaluating distance correlations. Besides the dipolar 
repulsion at short distances, an attractive force has been 
found at larger distances, and a minimum in the potential 
energy has been shown to occur at about 5.7 /an [7]. As 
a possible explanation, these authors propose that the 
deformation of the interface by the particles gives rise 
to a capillary interaction that varies logarithmically with 
their distance [7]. 

Quite generally, a capillary or elastic interaction is ob- 
tained when equilibrating an external force acting on the 
particles with the surface tension or the bending rigid- 
ity. A variety of such models have been studied, both 
for spherical and anisotropic defects [13-19]; these works 
deal with forces that act on the particle only, such as 
gravity. A somewhat different pattern arises for forces 



that operate both on the particle and on the surround- 
ing interface. For example, a polymer grafted on a mem- 
brane exerts a force at the point of attachment and an 
opposite entropic pressure on the interface [20] . 

In this Letter we study the deformation of a polar in- 
terface due to a charged colloidal particle, and we derive 
the resulting capillary interaction. The present work is 
confined to the case most relevant for micron size col- 
loidal particles, where the distance is much larger than 
the Debye length. After a reminder of the free energy of a 
deformation field, we calculate the electrostatic pressure 
profile exerted by a charge and the associated counteri- 
ons (as shown schematically in Fig. 1.) The electric-field 
induced capillary interaction is compared to experimen- 
tal findings and to the case of 5-force that was proposed 
in [7]. 

FIG. 1. a) Pressure profile n(r) induced on the interface 
by a charged particle (filled grey circle); open circles are the 
counterions. b) The corresponding deformation field u(r). 

We consider particles of charge Q trapped at an 
electrolyte-insulator interface. In the absence of forces, 
the interface S is flat and its energy reads 7 S, where 
7 is the surface tension. The charged particles and their 
counterions exert on the interface a pressure that is given 
by the normal component of the stress tensor and com- 
prises entropic and electrostrictive contributions, 



?r(r) =fc B T(n+ + n_)+7j 
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The first term is the entropic pressure of the ex- 
cess density of positive and negative ions, n± = 



,^e<t>(r)/k B T 



I), induced by N surface charges. 



Electroneutrality requires e f dV(n+ — n_) = 
remaining terms involve the Maxwell tensor 



EjEj 
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evaluated at the electrolyte and insulating side of the in- 
terface. Note that both the dielectric constant e and the 
normal component E z of the electric field vector are dis- 
continuous across the interface. The corresponding jump 
of the normal component of the Mawell tensor describes 
an electrostrictive force that arises from the difference in 
the electric field density in the two media. 

As a consequence of the inhomogeneous pressure, the 
interface is no longer flat; the deformation it(r) increases 
the total surface and thus the surface energy by the 
amount 
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where V is the 2D gradient in the interface plane. In the 
following, r denotes the in-plane coordinates, and z that 
normal to the interface. In the case of weak deformations, 
Vw| <C 1, we may expand the square root and retain the 
leading term only. Then the free energy functional reads 



/[«(')] = \\ dS(Vu) 2 - J dSn(r)u(r). 



(3) 



The actual deformation is determined by minimizing this 
functional. Linearizing the deformation about the equi- 
librium value u(r) and integrating the first term by parts, 
one finds that the mimimum free energy / occurs for the 
deformation satisfying the Young-Laplace equation, 



7V 2 w(r) + 7r(r) = 0, 



and reads 
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-\ J dS(Vu) 2 . 



(4) 



(5) 



As the most severe approximation of the present pa- 
per, we assume that the total pressure 7r(r) is the linear 
superposition 



7r( r ) = X^( r ~ r ° 



(6) 



where 7r(r) is the pressure profile of a single charge at 
the origin. In Debye-Hiickel approximation, the electro- 
static potential is a linear superposition of single-particle 
terms <f> — ^2 a f a , an d the total pressure is quadratic 
in the potential and its derivatives, ir(r) ~ <f> 2 . Besides 
the diagonal terms n a ~ (p 2 ^, it comprises cross-terms 
(p a (p/3 that are significantly smaller and thus have been 
neglected in (6). 



FIG. 2. The capillary interaction arises from the superpo- 
sition of the deformation fields of nearby particles. 

Because of the linearity of (4) , a similar relation holds 
true for the deformation 



w(r) 



u(r- 



It is then sufficient to consider two particles at positions 
±ip. Discarding terms that do not depend on their dis- 
tance p, one finds immcdiatly the pair potential 



A/(p) = - 7 j d 2 rVu (r + ^p) ■ Vu (r - ^p) . 



(7) 



The single-charge pressure 7r(r) comprises the entropic 
pressure of the excess soluted ion density and the electric 
field energy. It is obvious from Fig. I that the force on 
the charge is directed towards the electrolyte, whereas 



the electrostrictive force on the surrounding interface is 
in the opposite direction. It is essential to note that the 
total force vanishes, i.e., 
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dSn (r) = 0. 



(8) 



In view of (5) we need to calculate the gradient of the de- 
formation field. This is achieved most easily by applying 
Gauss' theorem on (4), 
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dsn ■ Vu 
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that relates the pressure integrated over the area D to 
the "flux" of the gradient field through the boundary 
3D. For a disk centered at the origin and polar coordi- 
nates r and 9, we have ds = rdO and n — r/r. Using the 
condition (8) and the fact that it is isotropic, the integral 
over the disk is replaced by that over the infinite space 
outside, 



Vw(r) 
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(9) 



Thus we have expressed the interaction free energy Af(p) 
by integrals over the single-charge electrostatic pressure 
7r(r) that remains to be determined. 

Starting from the screened potential ip(r, z) of a single 
interfacial charge at the origin, we evaluate the pressure 
profile according to Eqs. (I) and (2). Expanding the ion 
densities to quadratic order, n + + n_ = n s (eif/hsT) 2 , 
and noting k 2 = (2n s e 2 /eeifcsT), we find the entropic 
pressure 

k B T(n + +n_) = ^-K 2 <p(r,z) 2 . 

The relevant component T zz — l-(E 2 —E 2 ) of the Maxwell 
tensor depends on both normal and in-plane derivatives 
of the potential. The latter one, E r — d r ip, being contin- 
uous at the interface, the electrostrictive force simplifies 
to 



zz zz 



-El - £ -^E 



I ms TP; 



In Debye-Hiickel approximation (DHA), the 2D 
Fourier transform of tp is known, and explicit forms of 
the potential at the interface can be given in the do- 
mains separated by the Bjerrum length £b and the Debye 
screening length n^ 1 . Here we confine ourselves to the 
case of large distances well beyond k -1 and the particle 
size a, where the potential in the electrolyte reads [9-12] 



<p e i(r,z) = (Qe ins /2irK e el )e 



(r > «f 



(High charge densities may result in a renormalized value 
for Q.) The in-plane electric field E r ~ r~ 4 is neg- 
ligible at large r; for the normal component one finds 
E z \d = d z (p e \ = —Kifi e i. Its value at the insulating side 
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of the interface is obtained from the continuity condition 
£insE z \ ins = e ci E z \ ei . Thus the first term in T™ s - Tf z 
cancels the entropic pressure, the second one is negligible, 
and the third one gives 



7r(r) 



Eel 



Kip e i(r, 0) 



(10) 



Now it is straightforward to calculate the gradient of the 
interface deformation, 



V«(r) = 



1 



32tt 2 7 4k 2 r 6 ' 



(11) 



This remarkable result, Vu <~ r -5 , relies only on the 
asymptotic form of the potential ip c \ ~ r~ 3 and the fact 
that the net force (8) on the interface vanishes. As a 
consequence, the deformation field is negative and varies 
as«~ (— 1/r 4 ). 

Now we discuss the behavior of w(r) and Vu(r) at small 
distances. The singularities for r — > are clearly unphys- 
ical and would cause the deformation free energy A/ to 
diverge. The above screened potential ip e \ and electro- 
static pressure tt are valid at large distances r 3> a. 
At intermediate distances a < r < n" 1 screening is irrel- 
evant; evaluating the stress tensor with the bare electro- 
static potential ip c \ ~ r _1 , one finds a modified power law 
7r(r) ~ r~ 4 [21]. For micrometer size charged colloids, 
however, the particle radius a provides the most relevant 
physical cutoff. For this reason, it will be referred to in 
the remainder of this paper, although one should keep 
in mind that the actual situation may be more complex, 
especially for small molecules. 

Thus the power laws of both 7r(r) and Vu(r) cease to 
be valid at distances smaller than a. Moreover, Eq. (8) 
requires a strong negative force operating on the parti- 
cle at r = 0. An important conclusion can be drawn 
directly from the relation (8): Vu is positive everywhere 
and tends towards zero at r = 0. 

In order to regularize the surface integral in Eq. (7), 
we explicitly introduce a cut-off a and replace (11) with 



Vw(r) 
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32tt 2 1 e 2 cX n 2 (r 2 + a 2 ) 3 



(12) 



Although the precise form of this cut-off function is some- 
what arbitrary, it satisfies the limit Vu = for r — > as 
imposed by Eq. (9). For r 3> a, it shows the long-range 
behavior that has been obtained rigorously from (8) and 
(10). 

Inserting Vu(r±^p) in (7), one finds that the capillary 
interaction 

A/(p) = -( £ins g 2 /327r 2 7£c V 2 ) 2 / 
depends on particle distance and size through the integral 



/ drr / (19- 
Jo Jo i 



r 2 - \ P 2 



((r 2 + \p 2 + a 2 ^ 2 — r 2 p 2 cos 2 O^j 



This definite integral can be performed analytically. 
Since in all applications, the particle distance p exceeds 
significantly the size a, we retain only the leading term 
in powers of (a/p) and thus have 



2tt 



(l + 0(a 2 /p 2 )) 



The variation with the particle radius ~ a~ 2 does not 
depend on the precise form of the cut-off function chosen 
above. Thus we obtain the capillary interaction 



1 e 2 O 4 
Af(p)= ' £ ^ 



2 9 tt 3 e 4 «; 4 7 o 2 p 6 



(13) 



that is repulsive and varies with the inverse sixth power 
of the distance. Note that A/ depends on the particle 
charge, its size, and the Debye length. 

We compare our result for the capillary interaction 
with experimental findings and previous work. In Ref. 
[7], a minimum in the interaction potential of PMMA 
particles was reported to occur at a distance of 5.7 mi- 
crons. As a possible explanation, these authors consid- 
ered the competing electrostatic and capillary interac- 
tions. Intcrfacial charges are subject to the well-known 
repulsive potential 



V = 



gins V 



where p is dipolc moment formed by the charge and its 
screening cloud [10]. On the other hand, assuming an 
electric-field induced (5-force 7r(r) = — 7To#(r), one obtains 
an attractive capillary interaction Afs(p) = (ttq/2tt) In/?. 
The sum of V and A/,5 shows a minimum at finite dis- 
tance that has been considered in Ref. [7] . Note, however, 
that this force does not satisfy the condition (8). 

As shown schematically in Fig. 1 and discussed be- 
low (11), the electrostatic pressure profile induced by a 
screened surface charge comprises a long-range contribu- 
tion and a (5-force, 

7r(r) = g(r)r~ 6 - n 5(r), 

with a cut-off function g that is constant for r ^> a, van- 
ishes at the origin, and satisfies Eq. (8). By now it should 
be clear that the long-range contribution to the pressure 
completely changes the capillary forces: While that of a 
pure 5-force is attractive and depends logarithmically on 
distance, the capillary interaction of interfacial charges is 
repulsive and varies as A/(p) <~ p~ 6 . We conclude that 
the mimimum in the interaction potential reported in [7] 
cannot arise from V + A/. 

Finally we compare the relative magnitude of the forces 
arising from V and A/. If the particle's charge Q is 
located at the interface, the dipole moment reads p = 
Q/k, and the capillary interaction can be expressed in 
terms of the dipolar interaction V and the surface energy 
70 2 , 
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The total force — d p (V + A/) is always repulsive. A 
change in the exponent occurs where d p V = d p Af, defin- 
ing the cross-over distance 

._!{ Sins P 2 \ 1/3 

9 4U 2 47aV ' 

for p < p* , the capillary repulsion —d p Af ~ p~ 7 dom- 
inates, whereas at larger distances the dipolar force 
—d p V ~ p~ 4 takes over. For micron size particles at 
an oil-water interface with typical parameters, one finds 
that the cross-over distance hardly attains the microme- 
ter range and thus is of little relevance. Quite a different 
situation arises if a significant fraction of the charge is lo- 
cated at the side of the particle immersed in the insulator. 
For micron size particles, the electrostatic potential tp is 
enhanced by a factor of the order (na) 2 [6], and the capil- 
lary interaction increases by a factor («a) 4 . The modified 
dipolc moment reads p ~ Qa instead of Q/k, resulting 
in a larger cross-over distance p*. Note, however, that 
at distances much shorter than p*, the gradient of the 
deformation |Vu| is no longer small, requiring to go be- 
yond the quadratic approximation for the surface energy 
in (3). 

In summary, we have studied the capillary interaction 
induced by the electric field of charged particles at an 
electrolyte-insulator interface. The electrostatic pressure 
of a single charge consists of a force acting on the particle 
and the opposite electrostrictive force on the interface 
that varies as p~ 6 with the lateral distance. The resulting 
deformation, u ~ p~ 4 , induces a capillary interaction 
that is repulsive and obeys a power law A/ ~ p~ e . As a 
consequence, this capillary force cannot be at the origin 
of the secondary minimum observed recently for charged 
PMMA particles at an oil- water interface [7]. 

Note added. After submission of this Letter, a com- 
ment by Megens and Aizenberg on the paper by Niko- 
laides et al. [7] appeared in Nature [22]. In agreement 
with our discussion and our Eq. (8), Megens and Aizen- 
berg point out that there is no net force on the interface 
and thus invalidate the attractive logarithmic potential 
proposed in Ref. [7]. 
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